Interplay between ferromagnetism and superconductivity in tunneling currents 
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We study tunneling currents in a model consisting of two non-unitary ferromagnetic spin-triplet 
superconductors separated by a thin insulating layer. We find a novel interplay between ferro- 
magnetism and superconductivity, manifested in the Josephson efi'ect. This offers the possibility 
of tuning dissipationless currents of charge and spin in a well-defined manner by adjusting the 
magnetization direction on either side of the junction. 
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The coexistence of ferromagnetism (FM) and super- 
conductivity (SC) has recently been experimentally con- 
firmed [ij, I3 • This offers the possibility of observing new 
and interesting physical effects in transport of spin and 
charge. The spin-singlet character of Cooper pairs in 
conventional superconductors suggests that FM and SC 
are mutually excluding properties for a material. In- 
deed, coexistence of magnetic order with singlet SC is 
not possible for uniform order parameters Q. On the 
other hand, spin-triplet Cooper pairs 0, 1^ are in prin- 
ciple perfectly compatible with ferromagnetism. For in- 
stance, odd-in-frequency 3^=0 spin-triplet supercon- 
ductivity in superconductor-ferromagnet structures has 
been much studied in the literature The synthesis 
of superconductors exhibiting ferromagnetism, with si- 
multaneously broken U(l) and 0(3) symmetries, are of 
considerable interest from a fundamental physics point of 
view, and moreover opens up a vista to a plethora of novel 
applications. This has been the subject of theoretical re- 
search in e.g. Rcfs. 7], and has a broad range of possible 
applications. Also, focus on hybrid systems of fcrromag- 
nets and superconductors has arisen from the aspiration 
of utilizing the spin of the electron as a binary variable in 
device applications. This has led to spin current induced 
magnetization switching ^S^, and suggestions have been 
made for devices such as spin-torque transistors and 
spin-batteries |lflj |. Moreover, spin supercurrents have a 
long tradition in '^He 'll^, while recent work has focused 
on dissipationless spin-currents in unitary spin-triplet su- 
perconductors [1^ . 

This Letter addresses the case of two p-wave super- 
conductors arising out of a ferromagnetic metallic state, 
separated by a tunneling junction. Such states have 
been suggested to exist on experimental grounds, in 
compounds such as RuSr2GdCu2 08 UGe2 0, and 
URhGe and have been studied theoretically in e.g. 
Refs. [UlilEl- Coexisting FM and spin-triplet SC 
have also been proposed to arise out of half-metallic fer- 
romagnetic materials such as Cr02, and the alloys UNiSn 
and NiMnSb 0. We compute the Josephson contri- 
bution to the tunneling currents, both in the charge- 
and spin-channel, within linear response theory using the 



Kubo formula. Our assumption is that the superconduct- 
ing order is that of spin-triplet pairing, and we consider 
the analog of the A2-phase in '^He, i.e. SC order param- 
eters that satisfy |Ak||| 7^ l^kxxl 7^ and Akfj. = 0. 
In terms of the dk- vector formalism ^3] , we then have a 
non-unitary state since the average spin (Sk) = idk x dj^ 

^(|Ak||P — |AkxxP)z of the Cooper pairs is nonzero. 
Such a scenario is compatible with uniform FM and SC 
since the electrons responsible for ferromagnetism below 
the Curie temperature Tm condense into Cooper pairs 
with magnetic moments aligned with the magnetization 
below the critical temperature Tc. The choice of such a 
non-unitary state is motivated by the fact that there is 
strong reason to believe that the correct pairing symme- 
tries in the ferromagnetic superconductors (FMSC) dis- 
covered so far are non-unitary The exchange 
field will also give rise to a Zeeman-splitting between the 
t, i conduction bands, thus suppressing the SC order pa- 
rameter Ak|x 05 as illustrated in Fig. ^p). 

Another important issue to address is whether the SC 
and FM order parameters coexist uniformly, or if they are 
phase-separated. One possibility is that a spontaneously 
formed vortex lattice due to the internal magnetization 
m is realized in a spin-triplet FMSC However, there 
have also been reports of uniform superconducting phases 
in spin-triplet FMSC [IJ- A key variable determining 
whether a vortex lattice appears or not is the strength 
of the internal magnetization m 22]. Current experi- 
mental data on URhGe apparently do not settle this is- 
sue unambiguously, while uniform coexistence of FM and 
SC appear to have been experimentally verified in UGe2 
[23I 1. Furthermore, a bulk Meissner state in the FMSC 




FIG. 1: a) Tunneling of Cooper pairs between two non-unitary 
FMSC with non-coUinear magnetization, b) Band-splitting 
for I, J, electrons in the presence of a magnetization in z- 
direction, leading to a suppression of interband-pairing. 



2 



RuSr2GdCu2 08 has been reported in Ref. indicat- 
ing the existence of uniform FM and SC as a bulk effect. 
Consequently, we will use bulk values for the order pa- 
rameters and assume that they coexist uniformly. We 
emphasize that one in general should take into account 
the possible suppression of the SC order parameter in 
the vicinity of the tunneling interface due to the forma- 
tion of midgap surface states [2^ which occur for certain 
orientations of the SC gap. The pair-breaking effect of 
these states in unconventional superconductors has been 
studied in e.g. Ref. [2^. A sizeable formation of such 
states would suppress the Josephson current, although it 
is nonvanishing in the general case. Also, we use bulk 
uniform magnetic order parameters, as in Ref. [2^ . The 
latter is justified on the grounds that a ferromagnet with 
a planar order parameter is mathematically isomorphic 
to an s-wave superconductor, where the use of bulk val- 
ues for the order parameter right up to the interface is 
a good approximation due to the lack of midgap surface 
states. Moreover, we consider thin film FMSC such that 
the Lorentz-force acting on the electrons will be unable 
to accelerate particles in a direction parallel to the junc- 
tion. Our model is illustrated in Fig. ^). 

The main result of this Letter is that the Josephson 
current in the spin- and charge-sector between two non- 
unitary FMSC can be controlled by adjusting the rela- 
tive magnetization orientation on each side of the junc- 
tion provided that spin-triplet Cooper pairs are present. 
Our system consists of two FMSC separated by an in- 
sulating layer such that the total Hamiltonian can be 
written as i? = Hi, -\- + Ht, where L and R rep- 
resents the individual FMSC on each side of the tun- 
neling junction, and Hj^ describes tunneling of particles 
through the insulating layer separating the two pieces of 
bulk material. The FMSC Hamiltonian is given by [l^ 
i?FMSC = Ho + EkV'k-^kV'k, where Hq = JN-f {0)m^ + 
\ Ekff £k<T + Eka/3 ^ka/3^ka/3- Here, k is the electron 
momentum, Eko- = Ek — c^z, cr =T,J.= ±1, J is a spin 
coupling constant, 7(0) is the number of nearest lat- 
tice neighbors, m — {mx,my,mz} is the magnetiza- 
tion vector, while Aka/3 is the superconducting order 
parameter and foka^a = (c-k/jCka) is the two-particle 
operator expectation value. The ferromagnetic order 
parameter is defined by C = 2Jj{0){mx — ^my) and 
C2 — 2Jj{0)mz. We express the Hamiltonian in the basis 

V-'k = (ckT Cki cLk| cLkj^)''", where Cko- {cl„) are annihila- 
tion (creation) fermion operators. Note that there is no 
spin-orbit coupling in our model, i.e. inversion symmetry 
is not broken. Consider now the 4x4 matrix 

2Uidk-Ta,)t (ea-T-OV ^' 

which is valid for a FMSC with arbitrary magnetization. 
As explained in the introduction, we will study in detail 
a non-unitary equal-spin pairing (ESP) FMSC as illus- 
trated in Fig. QJl), i.e. Akfx = ^kll = 0, C = in Eq. 



(^. Since the quantization axes of the two FMSC are 
not aligned, one needs to include the Wigner d-function 
[igf denoted by 'D„'„{'d) with j ~ 1/2 to account for 
the fact that a f spin on one side of the junction is 
not the same as a t spin on the other side of the junc- 
tion. The spin quantization axes are taken along the 
direction of the magnetization on each side, so that the 
angle d is defined by niR • ihl = mnmi, cos{'&) where 
nii — |mi|. The tunneling Hamiltonian then reads Ht = 

Ekp..' ^-'-('?)(%4.'^p.' +^^kp4.'CkJ, where we ne- 
gleet the possibility of spin-flips in the tunneling pro- 
cess. The validity of the tunneling Hamiltonian approach 
requires that the applied voltage across the junction is 
small. Here, we will be concerned with the case of zero 
bias voltage, so that the tunneling Hamiltonian approach 
is appropriate. Note that we distinguish between fermion 
operators on each side of the junction corresponding to 
Ckcr and dka ■ Furthermore, we write the superconducting 
order parameters as Ak<TCT = \A]f:aa\e^^^^~^^'"'\ where R 
(L) denotes the bulk superconducting phase on the right 
(left) side of the junction while 6'k is an internal phase 
factor originating from the specific form of the gap in 
k-space that ensures odd symmetry under inversion of 
momentum, i.e. 9k = ^-k + t^- 

For our system, the Hamiltonian takes the form 
ffpMsc = ^^0 + Ha, Ha = J2ka ^L^ko-^ko-, where we 
have chosen a convenient basis cj)]^^ = (cj^.^, C-ktr) that 
block-diagonalizes ^k, and defined ^ktr = \\^£'ka^z + 

t^^Loaipx + iCTy) -)- l^\^„„(px - iCTy)] with Pauh matrices 
Oi. This Hamiltonian is diagonalized by a 2 x 2 spin 
generalized unitary matrix [/km so that the supercon- 
ducting sector is expressed in the diagonal basis ^j^.^ = 
^L^^ka = (7L'7_ka)- Thus, Ha =^Ek^ ^L^ka<^k^, 
in which Akcr = V-^aA-^aU^l = diag(£;kcr, --Ekcr)/2, and 

£^k. = v/eL + l^kaaP. 

In order to find the spin and charge currents 
over the junction, consider first the generalized num- 
ber operator Nap = Ek'^ka'^k/J- The trans- 
port operator in the interaction picture then reads 

P,,(i?)T*p4,(t)Ck^(i)e'*"^], where eV = fiL - MR is an 
externally applied potential. The general current across 
the junction can be written 

m = Y.Tap{Nap{t)), T = (-el, or) (2) 

a/3 

such that the charge-current is /*~^(t) = Io{t) while the 
spin-current reads I^(t) = (Ii (i), J2(i), /3(i)). Note that 
Eq. 10) contains both the single-particle (sp) and two- 
particle (tp) contribution. The concept of a spin-current 
in this context refers to the rate at which the spin- vector 
S on one side of the junction changes as a result 0/ 
tunneling across the junction, i.e. S = i[i?T,S]. As there 
is no spin-orbit coupling in our system, this definition 
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of the spin-current serves well [30j . The spatial compo- 
nents of I® are defined with respect to the corresponding 
quantization axis. We compute the currents by the 
Kubo formula, (iV„^(<)) = -i jl^dt' {[N^p{t), HtH')]), 
where the right hand side is the statistical expectation 
value in the unperturbed quantum state, i.e. when 
the two subsystems are not coupled. We now focus on 
the two-particle charge-current and z-component of the 
spin-current, such that only a — f3 contributes in Eq. 

Using linear response theory with the Mat- 
subara formalism, one arrives at {Naa{t))tp — 
2E<T^m{*a<T(eX^)e-^'*'=^}, where ^^aaleV) is 
obtained by performing analytical continuation 
iLOn — > eV + iO+ {uJn = 2'!mk^T,n — 1,2,3...) on 
^^„{iujn) = -/o^''^^dTe-"-(T,M„,(T)Af„,(0)), where 
M^,{t) = Ekp^-"(i?)Tkp4„(t)dp,(t), while T is the 
temperature, and T,- denotes the time-ordering operator. 
Explicitly, we find 

kp 

\.p=± 

(3) 

where E\^a = V^ko- + l^kaap, ^kcr = £ko- - Mr (R ^ L 
for k ^ p), and A^J^^^(ey) = lim,^„^ey+iO+ A[/(£^ka) - 
/(Ap£'pcr)]/[iwi/ -I- pE\^a - AE'pcr]; \p = ±1- 

In Eq. ||2J), we have used that T_k,-p = T^p which 
follows from time reversal symmetry, and f{x) is the 
Fermi distribution. Note that the chemical potential has 
been included in the excitation energies i?ka- In gen- 
eral, Eq. Q will give rise to a term proportional to 
cos(0^^ — (?2q)j the quasiparticle interference term, in ad- 
dition to sin(6'^g. — 9^^), identified as the Josephson cur- 
rent. In the following, we shall focus on the latter while 
a cornprehensive treatment of all terms will be given in 
Ref. ml. Consider now the case of zero externally ap- 
plied voltage {eV = 0). From Eq. (0), we see that the 
Josephson charge-current becomes 

= e E [1 + ^] co^(^p - MO];, - C) 

kpcra 

X |Tkp|'|Akaa||Ap,,|i^k7/(-Eka-Bp,), (4) 

with F^^ = E±[fi±Eua) - f{Ep,)]/{Eua T Ep^) while 
the expression for if ^ is equal except for a factor (— a/e) 
inside the summation. Observing that Eq. may be 
cast into the form Ij = Iq + Im cos{'d), we have thus 
found a Josephson current, for both spin and charge, 
that can be tuned in a well-defined manner by adjusting 
the relative orientation ■& of the magnetization vectors 
(For corresponding results in spin-singlet superconduc- 
tors with helimagnetic order, see Refs. 0112). Below, 
we discuss the detection of such an effect. 

In the limit where one of the superconducting order 
parameters vanishes internally on both sides, i.e. the 



equivalent of an Al-phase, we see that the interplay 
between ■(? and Oa-a- remains, as only one term contributes 
to the spin sum over {a, a}. In this case, the charge- 
and spin-current goes as cos^(z?/2) sin A^o-o-, where 
Ae^^ = 9];, - 9^, and Aq,, with a = {T, i}, q = {k, p} 
is the surviving order parameter. For coUincar magneti- 
zation = 0), an ordinary Josephson effect driven by 
the superconducting phase occurs. Interestingly, one is 
able to tune this current to zero for niL || — rriR (i? = tt). 

Another result that can be extracted from Eq. is a 
persistent spin- Josephson current even if the magnetiza- 
tions on each side of the junction are of equal magnitude 
and coUinear (i? = 0). This is quite different from 
the Josephson-like spin-current recently considered in 
ferromagnetic metal junctions There, a twist in the 
magnetization across the junction is required to drive 
the spin-Josephson effect. In this Letter, however, we 
have found a persistent spin-current in the two-particle 
channel even for coUinear magnetization. 

In the special case of = and equal SC phases on 
each side of the junction, i.e. 9];, — 9^,, Eq. reduces 
to the form /|^ = Josm'^{'&/2) sm{9^^ - 9f^) while 
/j^ = 0. This means that a two-particle spin-current 
without any charge-current can arise for non-coUinear 
magnetizations on each side of the junction in the 
absence of an externally applied voltage and with equal 
SC phases = 9^,; see, however, Ref. [s^ . 

It is well-known that for tunneling currents flowing 
in the presence of a magnetic field that is perpendicular 
to the tunneling direction, the resulting flux threading 
the junction leads to a Fraunhofer-like variation in the 
DC Josephson effect, given by a multiplicative factor 
sin(7r<I>/<I>o)/(7r$/$o) in the critical current. Here, 
— Trh/e is the elementary flux quantum, and $ is the 
total flux threading the junction due to a magnetic fleld. 
However, this is not an issue in the present case since 
the magnetization is assumed to be oriented according 
to Fig. QJl). Since the motion of the Cooper-pairs is also 
restricted by the thin-fllm structure, there is no orbital 
effect from such a magnetization. 

Note that the interplay between ferromagnetism and 
superconductivity is manifest in the charge- as well 
as spin-currents, the former being readily measurable. 
Since the critical Josephson currents depend on the 
relative magnetization orientation, one is able to tune 
these currents in a well-defined manner by varying i?. 
This can be done by applying an external magnetic field 
in the plane of the FMSC. In the presence of a rotating 
magnetic moment on either side of the junction, the 
Josephson currents will thus vary according to Eq. Q). 
Depending on the relative magnitudes of Jq and Im, the 
sign of the critical current may change. Note that such 
a variation of the magnetization vectors must take place 
in an adiabatic manner so that the systems can be con- 
sidered to be in, or near, equilibrium at all times. Our 
predictions can thus be verified by measuring the critical 
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current at eV = for different angles -d and compare the 
results with our theory. Recently, it has been reported 
that a spin-triplet supercurrent, induced by Josephson 
tunneling between two s-wave superconductors across 
a ferromagnetic metallic contact, can be controlled by 
varying the magnetization of the ferromagnetic contact 
[3^ . Moreover, detection of induced spin-currents are 
challenging, although recent studies suggest feasible 
methods of measuring such quantities j^^- Observation 
of macroscopic spin-currents in superconductors may 
also be possible via angle resolved photo-emission 
experiments with circularly polarized photon s 1361 . or in 
spin-resolved neutron scattering experiments |37l |. 

We briefly mention our results in the single-particle 
channel where we find that the charge- current and the 
z-component of the spin-current both vanish for eV = 0; 
see Ref. 0| for details. They are nonzero for eV / 
even if the magnetization vectors are collinear. We 
stress that the finding of a non-persistent z-component 
of the spin-current does not conflict with the results of 
Ref. ,2Zi, as their z-direction corresponds to a vector 
in the ccy-plane in our system. For Ako-o- 0, Isp(t) = 
2EkpEo^.2'^aW2?./3W|7kp|29m{o-0„A^'^i(-ey)}, 
and the component of the spin-current parallel to 
niL X mpi is seen to vanish for = {0, tt} at eV = in 
agreement with Ref. [ivj . 

We reemphasize that the above ideas should be exper- 
imentally realizable by e.g. utilizing various geometries 
in order to vary the demagnetization fields. One may 
also use exchange biasing to an anti-ferromagnet to 
achieve non-coUinearity ^38] . We have found an interplay 
between FM and SC in the Josephson channel for 
charge- and spin-currents when considering non-unitary 
spin-triplet ESP FMSC with coexisting and uniform 
ferromagnetic and superconducting order. 
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